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Abstract 

Consider a flock of birds that fly interacting between them. The 
interactions are modelled through a hierarchical system in which each 
bird, at each time step, adjusts its own velocity according to his past 
velocity and a weighted mean of the relative velocities of its superiors 
in the hierarchy. We consider the additional fact, that each of the 
birds can fail to see any of its superiors with certain probability, that 
can depend on the distances between them. For this model with ran- 
dom interactions we prove that the flocking phenomena, obtained for 
similar deterministic models, holds true. 



1 Introduction 

Consider a flock of k birds with positions and velocities denoted respectively 
by (xi, . . . , a;*;) := x and (v i, . . . , v^) := v. All individual positions and 
velocities are vectors in M 3 , and we refer to x and v as the position and 
velocity of the flock respectively. 

Cucker and Smale, in their seminal paper [I], propose to consider a dy- 
namic for this system in which each bird adjusts its velocity according to its 
own past velocity and a weighted mean of the relative velocities of the other 
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birds of the flock. More precisely, considering a discrete time framework with 
time step h > 0, at time t + h the velocity of bird i is determined by the 
velocities of all birds at time t by the equation 



Vi(t + h) = Vi(t) + a ij{ v j{ t ) ~ v i{t)) 

3=1 



(1.1) 



for certain weighting coefficients ay, i,j = l...,k. The dynamic of the 
system is completed by the innovations in the position, according to the 
equation 



To simplify the notation we write x[t] for x(th) and similarly for other func- 
tions of t. The interactions that Cucker and Smale [I] propose to consider 
are of the form 



for positive parameters K,a,(3, where |.| is the euclidean norm in M 3 . We 
remark that this proposal is related to Vicsek's model [10], and refer to [I] 
for a comparison between Vicsek and Cucker-Smale models. 

A modification of Cucker-Smale model was proposed by Shen [9], who 
considers a non-symmetric structure of interactions, where each bird is con- 
nected only with its superiors in a prescribed hierarchy. More precisely, Shen 
proposes in (9] to consider the same coefficients of the Cucker-Smale model if 
j is superior to i in the hierarchy and zero otherwise. This model was named 
Cucker-Smale model under hierarchical leadership. 

It is useful to consider that each bird in the flock constitutes a vertex of a 
graph. The interaction between two birds is modelled through the presence 
of an edge between the corresponding vertices. The hierarchical order imply 
that the adjacency matrix of the graph can be taken to be triangular inferior 
(see Definition [1] below). Our random selection of the interactions between 
birds imply that each entry of the matrix is a random variable. 

One relevant question for these models is to study the asymptotic be- 
haviour of the flock, more precisely, to answer under which conditions the 
birds do flock together, i.e. under which conditions all individual velocities 
converge to a common velocity and all relative positions of one bird with 
respect to another converge to a constant limit. We see that in this situa- 
tion, when birds flock, in a coordinate reference system that travels with the 



x(t + h) = x(t) + hv(t). 
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common limit velocity, each bird tends to a constant position, with asymp- 
totically vanishing velocity. 

Our proposal in the present paper, departing from the model proposed 
by Shen [9], is to consider that at each time step, the interaction between one 
bird and each one of its superior is random, including different possibilities: 
the presence and/or the magnitude of the connectivity can depend on the 
distance between them. The precise condition is f l3.ll) in Theorem [TJ See 
also Examples I4.lti4.3l In conclusion, we consider what can be named as the 
Cucker-Smale model under hierarchical leadership with random interactions, 
being our purpose to study the flocking phenomena for this model. 

Other works including randomness in the flock dynamic are [3J, where a 
random external perturbation is included at each step; and [5|, where the 
connectivity at each step and for each pair of birds has a constant, indepen- 
dent, failure rate (see example 14.11) . 

The rest of the paper includes Preliminaries in Section 2, our main result 
in Section 3, Examples in Section 4, the Proofs in Section 5, and Conclusions 
in Section 6. 

2 Preliminaries 

In this section we present our model and some general results for them. 

2.1 Flocks under Hierarchical Leadership (HL-flocks) 

Definition 1 (HL-flock, [9]). We say that a flock of k birds is under hier- 
archical leadership, or for short that it is an HL-flock, if the birds can be 
labelled as 1, . . . , k so that 

1. ciij (= a birdi led by birdj ) ^ implies that j < i, 

2. if we define the leader set of each bird i by 

JSf(z) = {j : a l3 + 0}, 
then for any i > 1, Jzf (i) ^ 0. 

Note that a flock is an HL-flock if and only if the adjacency matrix ((a^)) 
of the associated graph is lower triangular with at least one non vanishing 
entry at each row except the first one. Besides, bird 1 is the overall leader of 
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the flock and it moves with constant velocity. For i > 1, bird i can look only 
at birds with a smaller index, and should look at least at one of them. 

According to our definition, the velocity of bird i is given by the equation 

vt[t + 1] = Vi[t) + h M^M " V M)- (2- 1 ) 

In order to obtain a weighted mean of the velocities of the superiors of each 
bird in the previous formula, from now on, and through all the paper, we 
assume the following two restrictions on the parameters of the model: 

h < -—, < a l3 < 1. (2.2) 

The position of the flock is governed by (11.21) . The model considered in 
this paper assumes that ££{i) is fixed for each i = 1. . . . , k, and that the 
coefficients a^-, j G i = 1, . . . , k can be random. 

2.2 Change of coordinates 

In the models considered in [31111, the adjacency matrix A = (((%)) is sym- 
metric and in the proof of flocking a key role is played by the Fiedler number 
0, i.e. the second least eigenvalue of the Laplacian of the matrix A. Under 
hierarchical leadership, the model we consider, since the adjacency matrix 
A = ((dij)) is lower triangular, this approach is no longer useful. We refer to 
[9j for details. 

To overcome this difficulty we adapt methods from [HI [2] . Observe first 
that, as Jz?(l) = 0, it follows that vi[t] = «i[0] for all t, so, if we have 
flocking, all limiting velocities are equal to v i[0]. As a consequence of this 
fact, we decompose the velocities considering relative velocities with respect 
to vi[0]. For each bird % we write its velocity as 

Vi[t] = -«i[0]) + «i[0], 

and observe that the relevant dynamic for the properties we are interested 
in depends only on the differences Vi[t] := Vi[t] — «i[0], i — 1, . . . , k. Accord- 
ingly, the origin of our coordinate system will be the first bird, introducing 
then Xi[t] = Xi[t] — Xi[t], % = 1, . . . , k. 
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It is easy to see that the relative velocities Vi follow the same dynamic 
that the absolute velocities «j. In fact, from equation (12.11) we have 

Vi[t + 1] - «i[0] = Vi[t] -v 1 [0] + h a v MA ~ »<W) > 

= v i [t]-v 1 [0] + h a ij(Mt\-viM)-(v i [t]-v 1 [0})), 

and similarly for the positions. Therefore, we have the following dynamic 

Xi[t + 1] = Xi[t] + hvi[t], 
vt[t + 1] = +/jE i6 if(») <hj (vj[i\ - Vi[t]) . 

Despite the fact that V\[t] = and Xi[t] = for all t, sometimes we 
include them in the formulas for notational convenience. 

2.3 Basic results 

We now present some basic facts about the general model just described, 
in particular, in an HL-flock, the norm of the velocities, with an adequate 
norm, is non increasing in time. More precisely, for y = (yi, . . . ,yk) in M 3fe , 
denote 



IblL = max \\ Vi 

l<i<k 

Denote x = ||a;[0] II oo an d v o 



Lemma 1. Consider an HL-flock under the dynamic given in (12. ip under 
the constrains f l2.2p . Then: 

(a) The absolute velocity of the flock is decreasing in norm, i.e., for all 
t > 0, we have 

W + i]|loo<H*]IL- 

(b) The relative velocity is decreasing in norm, i.e., for all t > 0, we have 

H* + i]IL< Ht]\L. 

(c) The relative position satisfy 

\\x\t]\\oo < %o + hv t, 
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(d) The differences of the relative positions verify 

max WxAt] — xJtllloo < 2x + 2hv t. 

l<i,j<k 

The proof of (a) and (b) in this Lemma can be easily adapted from 
[2]. See also Dalmao [5]. The proof of (c) and (d) are a consequence of 
H*]||oc < ||u[0]||oo and (H2D. 

The norm ||x||oo can be considered as a measure of the dispersion of the 
flock and similarly for ||f ||oo- We then give the following definition. 

Definition 2. In a dynamical system governed by equations (12.31) . we have 
flocking when the following two conditions are fulfilled: 

(a) The relative velocities asymptotically vanishes, i.e. 

\\v[t]\\ao when t — > oo. 

(b) The relative positions converge, i. e. there exists x G M 3fc such that 

x[t] — > x when t — > oo. 

Remark 1. It is clear to see that a sufficient condition for flocking is 

oo 

53 Ht]|U < oo. (2.4) 

t=0 

In case of having a random system, a sufficient condition is 

J2mv[t]\\oo < oo, (2.5) 

t=0 

as this condition implies that condition (12.41) holds with probability one. 

3 Main Result 

A relevant issue in Cucker-Smale proposal is the mathematical formulation 
of the fact that, in order to prove flocking, certain control from below on the 
interactions between birds should be considered. With the introduction of 
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randomness in the model arises the question of how to mathematically for- 
mulate this fact in an intuitive way being also mathematically tractable. Our 
proposal in (13.11) is to express this condition through a conditional expecta- 
tion on the past of the system, where the vector T t = (#[0], v[0], . . . , x[t], v [t]) 
contains the information of postions and velocities of the flock up to time t. 

Theorem 1. Consider the HL- flock of Definition^ that follows the dynamic 
(12.31) . subject to condition (\2.2\ and 

^\t+mn)> (1 + M *_ xmr (3-D 

for j G S£[i\ % — 1, . . . , k, and for some p G (0, 1]. Then, if 

1) a < 1, or 

2) a = 1 and the coefficients 

^v^woFW e = 2 — k ' 

(where 1/0 = oo) satisfy the following conditions 

7 £ > k - i + 2, for£ = 3,...,k, and j 2 > k - 1; (3.2) 
we have flocking. 

Corollary 1. Under the conditions of TheoremUi for a < 1 and any deter- 
ministic sequence {5t} such that 

Vr¥- 2 exp (—, -4 zihtY-A < oo, (3.3) 

there exists a random time r such that 

Mt]\\oo < St, 

for t > r . 

Corollary 2. Consider the flock of Theorem^ Then, for a = 1, under the 
condition 

we have flocking. 
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4 Examples 



Example 4.1 (Random failure of connectivity). Consider that the presence 
of connectivity between each pair of birds is decided independently of the 
other birds and of what happened at previous times, with constant probabil- 
ity for all pair of birds and times. If the connectivity is present, the strength 
of this connectivity, i.e. the coefficient in the weighted mean in (11. ip obeys 
Cucker-Smale formula. The mathematical formulation of this situation con- 
sists in considering 

aiAt + 1] = ^ {l + WxM-xMW)* ' 

for some fixed a > 0, where the {X t i j} are independent Bernoulli random 
variables, with 

P(X flii = 1) = 1 - F(X tjij = 0)=p, 

for some p G (0, 1). It is clear that condition (13.11) holds. 

This same example is treated by other methods in [6]. A key role in the 
approach in [6] is played by the law of the iterated logarithm, having the 
advantage of giving a more accurate convergence rate. 

Example 4.2 (Random interactions). Assume that {X t i j} in (14. ip are non- 
negative, independent random variables that satisfy the condition 

infEX Mi >p. 

t,i,j 

It is clear that condition (13.11) holds. This is a generalization of example 14.11 

Example 4.3 (Connectivity in random environment). We select the active 
links for each bird at each time step considering that the probability of con- 
nection is a function of the distance between birds, and the strength of the 
interaction, in case of presence, is a constant. Hence transition probabilities 
of the system at time t + 1 depend on the vector of all the past values of the 
positions and velocities Tt- 

More precisely, for an HL- flock we consider coefficients t — 0, 1, . . . , 

with j G ^{i)i such that they can take the values 0, p with probabilities 

F( aij [t + 1] = p | T t ) = 1 - F( aij [t + 1] = 0\F t ) 

1 

- (i + \\ Xi [t] - Xj [i\\\r 

that verifies condition (13.11) . 
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5 Proofs 



We begin by an auxiliary result that resumes the contractive nature of the 
dynamic. Denote 

A = l + 2x , B = 2hv , w M = wm{\\v e [0]-v j [0]\\:jeSf(t)}, 
and assume that wqi > for all £ — 2, . . . , d (see Remark [2] below) . 
Lemma 2. Under the conditions of Theorem^ if a < 1 we have 



E 



t+i 



a=T+l 



< 



cxp 



hp 



1-aBo 



({Ao + Bot) 1 ^ - (A + Bot) 1 -*) 



and in case a — 1: 
E ' 



je:sf(/) 



<r=r+l 



< 



A + fe%T 
A + /MU M (t + 1) 



7* 



Proof. Based on properties of the conditional expectation, we have 



E 



t+i 



<J = T+1 



E 



E 



E 



cr=T+l 
t 



a £i [a] ) 



j'GJSf(€) (T=T + 1 

Assume that a < 1. Since JSf (£) 7^ and 

P 



E( aij [t + l]\F t ) > 



> 



• (5.1) 



P 



;i + ||x,[t]-x,[t]||)« " (A + BoO c 



we have 



E ( 1 ~ h Yl a tA t + 1 \ 



Ft) < 1 



hp 



which is a deterministic bound. Therefore 



(A + B ty 



E 



t+i 

n (i-h 



(T=T+1 



< 1 



< 



n 



(A + s t) c 

hp 

(A + B Q a)' 



E 



<T=T + 1 



iesf(<) 



exp 



< exp J — /ip 



(A + B a) 



(A + B a) a 



(5.2) 



Bounding the sum in the last line of (15.21) by an integral we obtain the desired 
bound. 

Assume now that a = 1. Analogously to the case a < 1 we have 



E 



aej[t + 1] 



> 



P 



A + h\\v e [0] -^[0]||t ; 



(5.3) 



therefore 
t+i 



E 



•Ft 



<r=r+l 



< exp 



P 



\-S)£ A ° +/ ' ii ^ [o] "^ [o]ik 



< exp 



-v 



Ho]-Vi[o] 



log 



^o + ^||^[o]--^[o]||(t + i; 

A + h\\v e [0] - Vj [0]\\t 
A + hw w r 



< 



A + hw i(t+l) 



This concludes the proof of the Lemma. 



• (5-4) 



□ 
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Proof of TheoremUl' Case a < 1 . We verify condition (12.51) and proceed by 
induction over £ = 2, . . . , k, considering the subflocks {1, 2, . . . , £}. 

Let us start with £ = 2. Recall that the first bird is the overall leader 
of the flock, in particular its (relative) velocity v\[t] = for all t. For the 
second bird, from Definition [1] we have Jz?(2) = {1}, then 

v 2 [t + 1] = v 2 [t] + ha 21 [t + l](vx[t] - v 2 [t]) 

= {l-ha 21 [t+l])v 2 [t\ 
t+i 

= Y[(l - ha 21 [r])v 2 [0}. 



(5.5) 



T=l 



By Lemma [2], we have 

E\\v 2 [t + 1]|| < Kexp 



(A + Bot) 1 - 



'1 - a m 



hp 



A 



l-a 



(l-a)Bo ^0 



this gives condition (12.51) . concluding 



where K = 1 1 [0] 1 1 exp 

the proof for £ = 2. 

We turn now to the general case. Assume that for j < £ there exist 
positive constants Cj such that for all t we have 



where k 



hp 



E\\vj[t]\\ < C^^e-^ ^ ^ 
From (12.31) . it follows that 



(5.6) 



(l-a)-Bo 

v e [t + i] = v t [t] + h a ^ + tfivAt] - MA) 

1-h y~] a ej [t + 1] \v t [t] + h Y a «[* + 

= (l-a[t + l])v t [t]+b[t + l]. (5.7) 
where we denote 

a[t + l]= Y a £i [t + 1] , b[t + l]=h ^ (5.8) 

Iterating (I5.7P and taking norms, it follows that 

t+i t+i t+i 

IW+l]||<II(l-Mr])|M0]|| + X;||6[r]|| J] ( 5 - 9 ) 



r=l 



T=l 



a=T+l 
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where we agree that n CT=t+2 (l — /ia[er]) = 1. By Lemma [2], there exists a 
constant C such that 



t+i 



E JJ(1 - ha[r)) < C e - K{Ao+Bot)1 



We have 



r=l 



E||6[r]||<fc Y, E(a^[r]||^[r-l]||) 
<h £ E||^[r 



and, by the induction hypothesis (I5.6p . for all t we have 

t+i 



t+i 



£e||6[t]||</i £ £E|K.[r-l]|| 

jej^(f) r=l 

t+1 

r i-2 p -K(A +Bo(r-l)) 1 - Q 



T = l 



je^) r=l 



t+1 



< C"(t + I)'" 3 e- K( - Ao+Bo ^- 1 ^ 



T=l 



where C = maxj^mCj. Using this information and Lemma [2j for the 
second term on the right hand side of (15. 9p . we have 



E 



t+i 



t+i 



£>[t]|| ; (l-ha[a\) 



T=l 



<T=T+1 



\l-3 



E 



t+i 



£llW 



T = l 



t+1 



' (i - M*]) 



<J = T+1 



t+1 



<C'(t+l r " 

T = l 



'^ e -K{(A +B t) 1 - a -(A +B T) 1 -«) e -K{A +B {r-l)) 1 



< K(t + lf- 2 e 



2-Kt 1 



for a convenient constant K. Then, the induction thesis (15.61) for I follows 
with C( = C + K. This concludes the proof of the Theorem. □ 
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Proof of CorollaryUl We have 

P(IH*]||oo > S t ) < tMvWU < ^C k t k - 2 e-^ Ao+Botl - a \ 

O t Of 

and, since the series in (13.31) converges, the result follows by the Borel-Cantelli 
Lemma. □ 

Proof of Theorem^ Case a — I. The proof is carried by induction as for 
the case a < 1. We begin by I = 2. 
From (15. 5p and Lemma [2] we have 

t+i 



E||v 2 [t + l]|| < |h[0]||EjJ(l-/iO2i[r]) 

< ||^2 [0 



T=l 

I 



An ^ 72 



A + hw oe (t+r 

In view of Remark [U under the condition 1 1 1^2 [0] 1 1 < p, the series of the 
velocities in ( 12. 5ft is convergent and we have flocking of the first two birds. 
Consider now the £-th bird. For the velocity, we have 

v t [t + l]=v e [t] + h a ^ + 1] i v AA ~ MA) 

= {1 - ha[t + 1]) v e [t] + b[t + 1], (5.10) 

where a[t + 1] and b[t + 1] are defined as in (15.81) . 
With this notation, we have 

t+i t+i t+i 

Mt + 1] = n ( x - + Yi n ( x - m*d, 

T=l r = l cr=r+l 

where we agree that II^h^I ~ ^ a [°"]) = 1- By Lemma [2] we obtain 

t+l . . , N 7£ 



E (1 - /mz[t]) < 

<X=T+1 



A + hw oe T 
A + /iu^ + 



Remark 2. T7ie case wne = is not considered in the previous computa- 
tions. It can be checked with slightly different arguments that in this case the 
convergence still holds. 
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Define recursively Si = A j£ — 1. In what respects b[t + 1] in (j5.10p . 
we have 

E\\b[r]\\=h ^Mr]\\vAr-l]\\)<h £ E||«i[r-l]|| 



An 



+ hw e (r - 1) 

<5«-i / 4 \ <5f-i 



< 



With this information, we conclude that 



E||^[t + 1]|| < ||v/ ft 



A + hw oe (t + 1) 







a + /iwoi(t + r 



where fc^ is a convenient constant that depends on the initial positions and 
velocities of the flock. If 8# > 1 we have flocking of the first I birds. It is easy 
to see that this condition is equivalent to (13. 2p with £ = k. This concludes 
the proof. □ 
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6 Conclusions 

The mathematical modelling of collective behaviour is relevant in many dif- 
ferent scientific frameworks. In our paper we study the existence of flocking, 
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i.e. the convergence to a common velocity and stable configuration in the 
positions of a system of k individuals that travel, adjusting their own veloci- 
ties as a wheighted average of the relative velocities of the other individuals, 
that, for simplicity, we name "birds". A very well known model in flocking 
was proposed by T. Vicsek, A. Czirok, E. Bet- Jacob, and O. Shochet in [TO] , 
where each bird is capable of seeing his neighbour birds within a certain 
fixed distance to compute the weighted average of velocities. Several simu- 
lation studies were carried out (see for instance [JJ) but the rigorous proofs 
of convergence results seems to be complex, as can be seen in the work by 
A. Jadbabaie, J. Lin, and A.S. Morse [7]. A fundamental advance in the 
modelling of the flocking phenomena was the introduction of power decaying 
(with the distance) coefficients for the computation of the weighted average 
by Cucker and Smale in [5]. This model allows to obtain conditions in terms 
of the initial state of the system in order to have flocking. A modification in 
the dynamic was proposed by Shen in (9] who considers a leadership system 
in the interactions, in which each bird has a set of superiors to obey when 
computing his future velocity (i.e. a subset of the flock), under the same 
Cucker-Smale power decaying interactions. It is clear that the interactions 
in Cucker-Smale model are stronger than the ones in the model proposed 
in [10]. Naturally arises then the question, whether it is possible to weaken 
the interactions of Cucker-Smale model in order to have more realistic mod- 
elling, preserving the flocking property. Our work inscribes in this direction, 
incorporating certain probability of failure for each connection at each time 
step and random interactions. Our main Theorem states that the flocking 
propoerty still holds under random interactions, provided that the condi- 
tional expectation of the interactions remains bounded from below by the 
Cucker-Smale interaction coefficient. It should be noticed that our proposal 
includes the possiblity of having at certain times a disconnected flock, but 
the probability of this event is handled in such a way that it does not prevent 
the existence (and mathematical treatment) of flocking. 
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